Transport properties of a molecule embedded in an Aharonov-Bohm interferometer 
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We theoretically investigate the transport properties of a molecule embedded in one arm of a meso- 
scopic Aharonov-Bohm interferometer. Due to the presence of phonons the molecule level position 
(sd) and the electron-electron interaction (U) undergo a polaronic shift which affects dramatically 
the electronic transport through the molecular junction. When the electron-phonon interaction is 
weak the linear conductance presents Fano-line shapes as long as the direct channel between the 
electrodes is opened. The observed Fano resonances in the linear conductance are originated from 
the interference between the spin Kondo state and the direct path. For strong enough electron- 
phonon interaction, the electron-electron interaction is renormalized towards negative values, i.e. 
becomes effectively attractive. This scenario favors fluctuations between the empty and doubly 
occupied charge states and therefore promotes a charge Kondo effect. However, the direct path 
between the contacts breaks the electron-hole symmetry which can efficiently suppress this charge 
Kondo effect. Nevertheless, we show that a proper tuning of the gate voltage is able to revive the 
Kondo resonance. Our results are obtained by using the Numerical Renormalization approximation 
to compute the electronic spectral function and the linear conductance. 



PACS numbers: 72.15.Qm, 72.25.Mk, 73.63.Kv 

I. INTRODUCTION 

Single-Molecular transistors are seen as basic elements 
to build nanometer-scale electronic devicesii"— Such 
nanostructurcs intermingle both the electronic and me- 
chanical degrees of freedom and therefore give rise to new 
physical scenariosi^"— Among their advantages, they of- 
fer a tremendous diversity and new potential function- 
alities. An archetype are the single-Cgo and -C140 tran- 
sistors in which single electron tunneling events are used 
to excite and probe the quantized mechanical degrees of 
freedom or vibrons of these moleculesj^i^ Recently these 
features have been also observed in suspended carbon 
nanotubesj ^^i^^ Importantly, these experiments shows a 
rather strong electron- vibron coupling which is consistent 
with the observation of the Franck-Condon blockade. 
The observation of a wide range of electron- vibron cou- 
plings in carbon nanotubes , i.e., from weak to strong 
couplings may motivate new experiments to test many 
of the theoretical predictions in vibron-assisted transport 
through molecules. Furthermore due to their flexibility, 
single molecules and especially carbon nanotubes can be 
contacted to non-metallic leads such as ferromagnetic or 
superconducting leads and they can even exhibit many- 
body physics such as the Kondo effect when they are 



strongly coupled to the reservoirs i^^'^ 

As we pointed out above the electronic motion through 
single molecules is sensitive to vibrational modes that af- 
fect dramatically the current Howi^^—^^^^ The center- 
of-mass motion mode modulates asymmetrically the cou- 
pling between the molecule and may open new channels 
for transport. The observation of multiple features in 
the current-voltage characteristic with identical spacing 
reflects the excitation of an integer number of a well 
defined vibronic mode4i^— Symmetric modes with an 
electron- vibron Holstein type of coupling may cause the 
Franck-Condon renormalization of the tunneling between 
the electrodes and the molecule. Such a mechanism sup- 
presses the sequential and cotunncling transport chan- 
nels through the moleculc^^. Importantly, the electron- 
vibron interaction leads to the polaron shift in which 
the charging energy U become effectively negative and 
the molecular level (e^) is renormalizedJ^i^ As pointed 
out in Ref. [l^ this intriguing aspect modifies the trans- 
port scenario opening a new transport channel through 
the molecular transistor. Owing to negative U, dou- 
bly occupied molecular states become energetically fa- 
vored. Thus, when the gate voltage is properly adjusted 
so that the energy of the double occupancy is tuned to 
the energy of the empty state, the molecule is at "reso- 
nance" for the tunneling of electrons in pairs^i^jl^ xhis 
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regime is known as the strong electron- vibron coupling 
regime. In chemistry, a variety of molecules in solution, 
have shown the feature known as potential inversion in 
analogy to the negative effective charging energy in the 
present context. In addition, the strong electron- vibron 
coupling regime can be achieved in vibrating carbon nan- 
otube quantum dots as some experiments have already 
probed. 

The functionality of such devices depends greatly on 
the ability to preserve phase coherence through the 
moleculej^"— Since the phase of electron wave function 
changes during the transport, it is of great interest to 
have access to the phase of the transmission amplitude 
in order to give a full characterization of the coherent 
transport through these devices. For this purpose we in- 
vestigate the widely known Aharonov-Bohm (AB) effect 
as a valuable tool to investigate quantum coherence of 
electrons. When the coherence of a circulating electron 
wave packet enclosing a magnetic flux $ is preserved, 
this results in an extra flux-dependent phase shift. In 
the simplest realization of an AB interferometer, an in- 
coming electronic wave function splits into two paths, 
which join again into the outgoing electronic wave func- 
tion. Applying a magnetic flux which threads this closed 
geometry, the outgoing wave function acquires a flux- 
dependent phase, (/? — e^/hc, where $ = BS is the flux 
{B is the applied magnetic field, S is the enclosing sur- 
face and e is the electron charge). As a consequence, the 
transmission is a periodic function of ip. The evolution 
of the electronic phase in AB geometries (ring and ring- 
like structures) was investigated in a series of pioneering 
works Their theoretical predictions were tested in 

many type of experimentsi^""— 

In our particular case, we consider a molecular state 
occupied by an unpaired electron inserted in one of the 
arms of the AB interferometer. The other charge states 
with lower energies are doubly occupied and are there- 
fore not active transport channels. The other arm acts 
as a direct path between the two electrodes [see Fig. [T] . 
This situation corresponds to the case of single-molecules 
with a non vanishing spin that behaves as magnetic 
impurities^ exhibiting many-body physics such as the 
Kondo effectJ^ At sufficiently low temperatures, the lo- 
calized spin residing in the molecule is able to form a 
many-body spin singlet state with the delocalized elec- 
trons in the electrodes. In this manner the molecular 
spin is effectively screened and a narrow quasi-particle 
resonance emerges in the local density of states (DOS) 
at the molecule site. The width of such a singular- 
ity is related with the natural energy scale of this phe- 
nomenon given by the Kondo temperature, Tx- Thus, 
Tk is a measure of the binding energy of the Kondo 
many-body singlet state. In transport through quan- 
tum dots, the Kondo effect leads to the unitary con- 
ductance as predicted theoretically^ and demonstrated 
experimentally4i The unitary limit of the conductance 
in the Kondo regime is related to a transmission phase S 
equal to 7r/24^ This prediction has been tested by em- 




FIG. 1. (Color online), (a) Schematic of the molecule embed- 
ded in an Aharonov-Bohm interferometer, (b) Energy level 
diagram for the molecule. 

bedding a quantum dot in one of the arm of an Aharonov- 
Bohm ring4ii^ Theoretical studies addressing this prob- 
lem have revealed the strong sensitivity of the phase shift 
to Kondo correlations. In our particular sys- 
tem we expect striking and interesting transport features 
due to the interplay of the electron-phonon (e-ph) cou- 
pling and Kondo correlations. In the Coulomb blockade 
regime, experiments performed with molecules such as 
Ceo and C140 showed an internal mode of vibration of 
huQ 5 — lOmeV— while the vibron mode is of order 
1 meV in carbon nanotubesJ^ This energy is compara- 
ble or even larger than the Kondo temperature observed 
in these molecules which is of the order of few Kelvin. 
Additionally, for some molecules the Coulomb charging 
energy becomes considerably reduced due to the effec- 
tive screening coming from the electrodes. In this way 
the electronic and some vibrational energy modes (typ- 
ically of the order of 0.01-0.1 eV) may become compa- 
rable. This is for example the case where the interplay 
between the Kondo effect and the vibrons leads to a rich 
physical phenomena reflected in the molecular transport 
as shown below. 

This paper is organized as follows, in Sec. |ll] we in- 
troduce the model Hamiltonian. We discuss different 
transport regimes depending on the strength of the e-ph 
coupling, the local mode frequency, luq and the electron- 
electron (e-e) on-site interaction U. In sec. IIIIl we briefly 
present our theoretical approach based on the numer- 
ical renormalization group (NRG). Expressions for the 
current and linear conductance are derived in Sec. IIVI 
Our numerical results are detailed in Sec. |Vl Finally, we 
conclude this work by summarizing our main findings in 
Sec. lVIl We add two Appendices that help understanding 
the work. In Appendix A, we obtain the effective Hamil- 
tonian for the different transport regimes. Appendix B is 
devoted to calculate explicitly the molecular Green func- 
tion and the phonon Green function using the equation- 
of-motion approach. 
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II. ANDERSON-LIKE AND KONDO-LIKE 
HAMILTONIANS 

In the following we introduce the model Hamiltonian 
in order to describe the molecular transport through the 
Aharonov-Bohm interferometer. In a first step, we dis- 
cuss the model we use to describe a isolated molecule. 
Next, we insert the molecule in the AB interferometer 
and analyze how the Kondo correlations emerge. 



A. Isolated molecule coupled to a single local 
phonon mode 

In order to describe the isolated molecule in the pres- 
ence of vibronS)2^ we use the following model Hamilto- 
nian 



n 



Hoi 



M 



T~Le 



-ph 



ph 



where 



a 

T-Le-ph = A(a + a)) [ud - 1) , 
T-Lph = LOQa)a. 



(2.1) 



(2.2a) 

(2.2b) 
(2.2c) 



The isolated molecule is described by the hamiltonian 
T-Lm'- The operator d\ (d^) is the creation (annihilation) 
operator of one electron in the unique active molecular 
orbital level (sd) for transport. U is the on-site Coulomb 
interaction, and Uda — d^d^ is the quantum occupation 
of the molecular level per spin [ud = '^^nda)- a^(a) 
is the creation (annihilation) operator for an elemen- 
tary quantized mechanical excitation or vibron which we 
view here as a local phonon mode with a frequency wq- 
In order to simplify our description of the e-ph interac- 
tion we only account for single mode for phonons cor- 
responding to any of the vibrational normal modes of 
the moleculeJ^"— This is the extensively used Holstein 
model to describe in a simple fashion the observed molec- 
ular energy spectrumi^^ This is a good approximation as 
long as the step-like features in the I-V characteristic of 
some molecules are attributed to excitations of a well de- 
fined single mode4i^— The generalization to more modes 
is straightforward. 



1. A first insight on the effect of the electron-phonon 
coupling 

In order to gain some first simple insight about the ef- 
fect of the vibrational modes on the electronic properties, 
it is convenient to describe the phonon degrees of free- 
dom by a simple harmonic oscillator potential coupled to 



the electrons. This can be written as 



V{x) = -mulx^ - Xla^ + a){nd - 1) 




V2\{7ld - 1) 



— A(n,-1)) -^{nd-l)\ 

(2.3) 



where 



1 



1 

\/2mwo 



[a) + a) 



(2.4a) 
(2.4b) 



Thus, the effects of the e-ph coupling leads to both a 
shift of the average position of the harmonic oscillator 
and an additional negative energy shift corresponding to 
the energy gain due to polaron-formation. By rewrit- 
ing this polaronic shift, we can anticipate the resulting 
renormalizations: 



Wo 



)2 y2 

{nd-lf = nd{nd-l)^ . (2.5) 

cjQ Wo 



Reabsorbing these terms into the original Hamiltonian, 
we obtain the renormalizations 



off 



u ^ 



2A^ 

Wo 



Ed, off 



>2 

-. (2.' 
Wo 



Note that the potential minimum position now depends 
on the occupation number of the dot. The renormal- 
ization of the e-e interaction gives rise to very differ- 
ent physical scenarios depending on the e-ph coupling 
A. From Eq. (12. 6p we see that Ucs can be negative 
when 2A^/(t/wo) > 1. This corresponds to the so-called 
strong e-ph coupling regime where the e-ph interaction 
induces an effective attraction between electrons. On 
the other hand, for a weak e-ph coupling, i.e., when 
2A^/(J7wo) < 1, the bipolaron quasiparticle is not en- 
ergetically favorable since the induced effective e-e inter- 
action remains positive. 



B. Including the leads: Kondo correlations 

Let us assume that a molecule is attached to the source 
and drain electrodes and placed in the lower arm of an 
Aharonov-Bohm interferometer (see Fig. [1} . The upper 
arm connects directly the source and drain contacts by a 
direct path. The simplest Hamiltonian able to describe 
the system is the Anderson-Holstein Hamiltonian: 



■H = Hc + HHol+nT, 



(2.7) 



4 



where 



Oi—L I R,k,(7 
a=L/R.,k,cr 



(2.8a) 
(2.8b) 
(2.8c) 



where we have defined i 



is the creation 



"L{R),k,cr FL(_R,),fc,o 

(annihilation) operator for an electron in the state k with 
spin a in the lead L{R). We have assumed in our analy- 
sis that the leads are non polarized such that the disper- 
sion relation is spin independent. Vq, is the tunneling 
amplitude for the a-electrode with the molecule. W is 
the amplitude probability for electron transfer from the 
left electrode to the right one, and — 27r$/(/)o is the 
Aharonov-Bohm phase, with 0o = he/ e being the quan- 
tum flux. 

Kondo physics arises when the molecule is strongly 
contacted to the reservoirs. "^^^ High-order tunneling pro- 
cesses between the delocalized electrons in the contacts 
and the localized electron in the molecule lead to screen 
effectively the molecular spin and to create the Kondo 
state. Then, a spin singlet many-body state is formed by 
the conduction band electrons and the localized unpaired 
electron in the molecule. As a result, the molecular DOS 
shows a very narrow peak at the Fermi energy [Ep)^ 
This is the so-called Kondo resonance. In order to see 
how Kondo physics emerges in the presence of e-ph in- 
teraction, let us first perform a series of transformations 
of -H. 

The first step consists in eliminating the e-ph coupling 
in T-Lm by applying the Lang-Firsov canonical transfor- 
mation^'^'^^ to the full Hamiltonian % — t- U^T-LU where 
the operator is defined as 



exp 



-{nd - l){a) - a) 



(2.9) 



The resulting Hamiltonian % — He + T~Lm +'Ht reads 
in the parity basis 



l—ej o.k.cr 



(2.10) 



k,k'{T 



k'a^oka I , 



Um = LOoa'^a + [ Ed + —] Ud + [U - — I nd-\nd\, ■ 



Ut = V2F E 



•■,k,a L 



-A(at-a)/<.o sin((p/2)4^d, + /i.e. 



where we have defined the even-odd parity basis as 

Cefc. = (e'^^/^CLfc. + e-^'^'^Rk.) , (2.11a) 

(2.11b) 



and Vl = Vu = V. Notice that now the vibronic degrees 
of freedom appear in a nontrivial manner in the tunneling 
part of the Hamiltonian. As we anticipated in Sec. Ill A 11 
the charging energy and the molecular level position are 
renormalized by the phonon frequency and the e-ph cou- 
pling as C/cff = 11-2X^/000 and £d,off = £d + A^/wq- 
Due to this renormalization Ucs can be negative induc- 
ing a bipolaronic attraction between the electrons (when 
2X'^/{uJoU) > 1). In this case and when A = Sd + U = 
(that corresponds to the particle- hole or symmetric case), 
the zero and the doubly occupied molecule states are de- 
generate and the low-energy excitations consists of charge 
fiuctuations with a large gap for the spin fluctuations. 
These energy excitations can be described in terms of 
an effective Kondo model in the charge sector where the 
molecular charge fluctuates between = and Ud = 2. 
On the contrary, in the weak e-ph coupling limit (when 
2X^/(10011) < 1) the magnitude of the e-e interaction is 
effectively reduced but it does not favor the bipolaron 
formation. Here, the physics corresponds to the common 
spin-1/2 Kondo effect where the singly occupied molecu- 
lar state is the ground state. Let us detail both regimes. 



1. The weak electron-phonon coupling limit: 
2\^/{ujoU) < 1 

To treat this limit, we first perform a second order 
Raileigh-Schrodinger perturbative calculation (or projec- 
tion method) in the tunneling Hamiltonian I-Lt, we derive 
an effective Kondo-type Hamiltonian of the form^ 



n 



Hoi + ^12-^; ^21 



£■10 — Eom 



,(2.12) 

The eigenenergies of the isolated molecule are given in 
Appendix A. The projectors read 

noi^V2 J2 Vi{m\U-\0)cl^J,il-nda), (2.13a) 
^21 = ^2 E ^/(™|W+|0)sgn(a)4ndsQfc,, (2.13b) 

where we have introduced 

= exp [±A(at - a)/a;o] , (2.14) 

and defined a — —a, Ve — cos{(p/2)V and Vo = 
sm{ip/2)V. The explicit calculation of Heff is given 
in Appendix A. After some algebra the effective Kondo 
Hamiltonian can be rewritten in a more standard form 



5 



as [see Eq. (jAll|) ] 



2. Strong electron-phonon coupling: 2}? /{ljjqU) > 1 



af3 



5" 



(=1 



I J2 K^^ci^^cp,,, (2.15) 



where the first two hnes are the Kondo interaction part 
while the last line corresponds to the potential scattering 
part characterized by the coupling K"'^ . Notice that the 
Kondo part is described by the isotropic Kondo model 
that is characterized by a single exchange coupling [see 
Eq. (jTzl) and Eq. lfMj) in Appendix A]: 



V*Vp\{m\U-m-' 



+ Ed - mujo 



A2 



^ Vo.Vp*\{m\U+\0)\^ 



(2.16) 



with 



{m\U^\Q) = exp 



.If A 

2 V^o 



A 

± — 
I V Wo 



|m) = -L=(at)"|0). (2.17) 



To shorten the notation we omit the sub-index ph in 
the molecular states. From Eq. p.l4p one has U~ = 
(— 1)™Z^+ and for the symmetric model, Sd = —U/2, 
Eq. p.iep becomes: 



2A^ I 2mujQ 



where 

r° ^ Jcos^{(p/2), J°° = Jsin2((p/2), J" 
with 

J 



(2.18) 



^ sin(v3) , 



8V^^_ |(m|Z^-|0)|2 



2A^ 1 2mu]o 



(2.19) 
(2.20) 



Expanding J"'' (A) around A = wc find 



j"^(A) « -V^Vp 



1 + 2 1 



U 

2tJo 



(2.21) 



Then, using the definition for the Kondo temperature in 
the isotropic Kondo model one has 



Tk oc exp 



1 



^(A)po 



(2.22) 



As a result, the e-ph coupling leads to an increase of the 
Kondo temperature in the weak e-ph regime as already 
found in Ref. (l^l by Cornaglia and coworkers. 



In this regime the e-ph interaction induces a bipo- 
laronic attraction between electrons since [/off becomes 
negative. In this case the ground state is a molecular 
doublet state composed of the doubly occupied molecu- 
lar state |2, 0) and the empty molecular state |0, 0) which 
are degenerate. The low-energy excitations thus consist 
of charge fluctuations and can be described by an ef- 
fective Kondo model in which the role of spin is played 
now by a pseudospin variable that represents the two 
states of the doublet. Importantly, there is no rotational 
invariance in pseudospin space and therefore the effec- 
tive model corresponds to an anisotropic Kondo model 
with different couplings, namely Jy and J^. In order to 
gain some insights about the physical properties of this 
regime, we employ again the projection method. Before 
this, it is convenient to map the negative- J7 Anderson- 
Holtstein model to an equivalent model with positive in- 
teraction. This can be carried out by using a particular 
particle- hole (p-h) transformatioui^SiSL First, we intro- 
duce new operators d, c via 



C-Lki- 



CRki 



Ikl 



(2.23) 



The p-h transformation yields Eda — > £da ~ U/2 + 
sgn(cr)A/2 and U ^ U = -U where A = 2ed + U and 
Ska = ~^ka- assuming the symmetry condition for the 
conduction band. The tunneling amplitudes are given by 



(2.24) 
(2.25) 



The p-h transformation maps the electronic states of the 
molecule to the ones of the original model in this manner 

10).^ Il), 12).^ It), U)^|0), |t)^|2). (2.26) 

Now we perform the polaronic transformation^^'^^ com- 
bined with a canonical transformation such that the full 
Hamiltonian in the parity basis [see Eq. (I2.1ip ] reads (see 
Appendix A) 

l—e/o,k,a 

+ Sgn((T)M^ (cL'o-Cefe^T - cl^'a^oka 

k,k' ,a 

/ ;\2 \ _ _ /_ 2A^\ 

TiM = Kda I d^da + C/ H I UdfUdl 

^ V ^oj V ^0 J 

+ cjoo^a , 

J2 (Ke-^s"(")^/'^°('^'-")cos((^/2)4^d<, 

f_—e/o.k.(T 

+ v^e-'sn{.)xMa^-a) sin((^/2)4^d, -I- /i.e.) , (2.27) 
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where we assumed V\iia = Kr- The final step consists 
in obtaining an effective Hamiltonian by using again the 
projection method. We employ Eq. (|2.12|) and we refer 
to Appendix A for details 



(2.28a) 

7^21=72 V,*fe.("^|Z^+'^^"^"^|0)sgn(a)dtHdsC,fe,, 

(2.28b) 



l.k.i 



with a = —(7, and Vea = cos((^/2)Vcr and Voa — 
sin((^/2)VJj. After the projection, the effective Hamilto- 
nian He/ / corresponds to the anisotropic Kondo model 
with the following exchange couplings (see Appendix A) 



J, 



Vc.Vp\{m\U-\Q)\'' 



1 



U 



/3 _ 8 V^Vp{0\U-\m){m\U-\0) 



2A2 



- 1 



U 



LUoU 

V^Vp\{m\U-m^ 



Ml. 



1 



U 



(2.29a) 
(2.29b) 
(2.29c) 



The Kondo temperature in this case is given by^ 

/ J, \ l/'/llPo 

Tk'xD y-j-J . (2.30) 

In the limit of large A/wq, J_l/ J|| behave as 

Jj./ J|| « exp [~2{\/ujof] . 
Then, the Kondo temperature Tk reads^ 



(2.31) 



Tk oc exp 



r 



(2.32) 



In contrast to the weak e-ph coupling case, the Kondo 
temperature decreases sharply with increasing A (F = 
with pq = pi^ = pf( being the DOS of the con- 
duction band electrons, we took the same for both con- 
tacts). Equations (|2.22p and (|2.32p are the main results 
of this section. 



NRG method is employed to solve the full Hamiltonian. 
This allows us to compute the linear transport properties 
through the AB interferometer such as the linear conduc- 
tance. It is known that NRG method suffers from loss 
of accuracy coming from the discretization of the con- 
duction band and the truncation of high energy states. 
This problem is serious in our system because the Hilbert 
space grows exponentially with the number of conduction 
channels in the leads. Any careless truncation can lead 
to significant deterioration of thr results. We thus exploit 
recent improvements in the NRG procedure as well as a 
symmetry property of the system to obtain high-quality 
predictions)^ At very low frequencies the NRG schemes 
developed by R. Zitko and T. Pruschke in Ref. [g^ and 
W. Hofstetter and G. Zarand in Ref. [65| yield similar 
results. These results are in agreement with those ob- 
tained by employing directly the Friedel-Langreth sum 
rule''^ valid only at zero temperature. However, for the 
high frequency spectrum it seems more convenient to em- 
ploy the NRG code, for instance if we were interested in 
the vibronic sidebands, developed in Ref. (6^. Since we 
are interested in the linear transport at zero temperature, 
then we simply use the Friedel sum rule that agrees with 
very refined NRG schemes in the low frequency regime. 

According to the standard NRG procedure, Eq. (|2.7p 
can be rewritten as (note that the conduction band is 
now discrete with states denoted by the /-operators): 

= Hm + T~Le-ph + T~Lph 



+ E V (cOs(i^/2)4/oea + S\w{Lp /2)d\fooa + h.c) 

(7 

+ E EE^^"^'c„(/L./«+i«^ + M 

a=e/on=0 

+ - V^E (-^Oecr/oeo- - /ooct/ooct) , (3.1) 
(J 

with ^ = tt'^W'^p^ and 

fnea = ^ (e^'^'^fuL. + C-'^^ fr^R^^ , (3.2a) 

= (-e''^/2/«L. + e-^'^/^nfl.) , (3.2b) 

the parity basis. Now we apply the NRG method to solve 
the Hamiltonian. We define a sequence of Hamiltonians 



III. THE NUMERICAL RENORMALIZATION 
GROUP APPROACH 



In order to describe both the spin Kondo effect, oc- 
curring for 2X'^ / {coqU) < 1, and the charge Kondo ef- 
fect, occurring for 2X^/{ujoU) > 1, on the same footing, 
we have therefore used the Numerical Renormalization 
group (NRG) techniquei^i^ Unlike other approaches, 
the NRG is not restricted for the values of the param- 
eters (A, Wo, U, F, W, (f, etc). 
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N- 



+ ^ J— (cos(^/2)4/oe. + sin((^/2)4/ooa + h.c.) 



a=e/o "=0 a 



is the background transmission with ^ = tt'^W'^Pq. The 
symbols 5 and 3? denote the imaginary and real part of a 
complex number and Gda da the quasi-particle retarded 
Green function for the localized electrons at the molecule 
site. The background reflection is obtained immediately 
by TZb — 1 — It is also useful to define F = r/(l + 
a = ATLTn/T^ with T = r^ + rfl (r^f^) = 7rpLiB.)Vl(^Ry 
below we discuss only results for the case of symmetric 
tunneling rates, i.e., Tl = = F = TrpoV^ with a — 1). 
In the linear response regime, Eq. (|4.2p can be written as 



E (fLh^^ - fLfoo.) ] , (3.3) ^ = 1^-0 ? E / * rAe)Ue) , 



(4.5) 



that follows the recursion relation 

Hn+I = VA-Hn + ^ E [fhaafN+laa + h.C. 
OL—e/o <y 



with 



(3.4) 



where /eg is the equilibrium Fermi function. At fcsT = 0, 

%{e)6{EF - 6) . (4.6) 



Thus, 



-Hn = 



1 + A- 



-n 



N ■ 



(3.5) 



From consecutive diagonalizations of the Hamiltonian, 
we obtain the eigenvalues and eigenstates. With them we 
build the molecular and phonon Green functions that are 
needed for the calculation of the linear transport prop- 
erties. In order to improve the quality of these Green 
functions, we follow Jeon et aLr- and relate the molecule 
and local phonon Green functions by using the equation 
of motion (EOM) technique. We relegate the details of 
such a calculation to the Appendix [Bl 



IV. LINEAR TRANSPORT: CONDUCTANCE 

The electrical current through the left barrier is calcu- 
lated through the simple relation 



lL^-eNL = —[H,NL], (4.1) 

where = J2kcr ^Lka^Lka- After some algebra the 
current through the AB interferometer takes the usual 
Landauer-Biittiker formi ^^i^^ 

^= JE/ deTAe)[fLie)~fH{e)], (4.2) 

where fL{R.) is the left (right) Fermi function and the 
transmission probability reads 



^% + v/«W6Cos((^)r5R [^5.,d.(e)] 



- {a [1 ~ Ttcos'i^)] - Tfc} F3 [g^,,^,M , (4.3) 



where. 



4^ 



(1 + 0= 



(4.4) 



G = ^^r.(i;F), atkBT = 0. (4.7) 



This can be further simplified by taking the quasi-particle 
Green function for the localized electrons written in a 
Dyson type equation as follows 



1 



e-£d-S(e) 



(4.8) 



At fcsT = 0, in the pure Kondo regime, the imaginary 
part of the interaction self-energy evaluated at the Fermi 
energy vanishes, then we have 

^[^{Ef)] = -^, q=l{ed + m{EF)]} ■ (4.9) 

Using this notation we can write a more transparent ex- 
pression for the transmission probability in a generalized 
Fano form as: 



T{Ef) = % 



with a Fano parameter given by 



cot(5,i 



oRb 

IT 



cos{ip) . 



(4.10) 



(4.11) 



Notice that we have introduced Snon as a non-resonant 
phase shift owing to scattering of electrons through the 
direct path. Alternatively, in terms of resonant and non- 
resonant phase shifts the transmission probability can be 
written as: 

T{Ef) ^ [Tb + aTZb COS^ {(p)] Sin^ {Sres - Snon) 

I • 2/ \ • 2/r \ T- sin {Sres ^ Snon) 

-I- asm (^)sm (dres) ^ Tb r^— r 

sm (d„on) 

-I- asin^(i^) sin^(^res) , (4.12) 
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FIG. 2. (Color online) Charge susceptibility Xc{^) ^ func- 
tion of the e-ph coupling X/uq. Parameters are Ed = 0.05, 
U = 0.1, 2r = 0.0016, ip^O,Tb = 0, keT = 0. 



where the resonant scattering phase shift Sres can be re- 
lated to the quantum average occupation number Ud — 
J2a(^td'cr} using the Friedel-Langreth sum rulei^ 



= tan 

2 



r 



(4.13) 



V. NUMERICAL RESULTS 

In the following we present our numerical results based 
on the NRG method. We fix = 0.1, 2T = 0.016, 
Wo = 0.05 and fc^T = and we vary the rest of pa- 
rameters. We have chosen a value of the phonon vibra- 
tional frequency ujq of the order of the Coulomb energy 
scale which corresponds with the experimental situation. 
Using these parameters the bipolaronic attraction occurs 
for X/luq > 1 whereas the weak e-ph coupling limit takes 
place for X/ujq < 1. 

We start investigating how are affected the Kondo cor- 
relations by the presence of phonons (A 7^ 0). In order to 
do so wc expand Eq. (j4.12l) around the symmetric point, 

= — C//2 (for % = and (p = 0), we have 



G{ed)^G{s*d)-Go{^^y{sd-s*,)' 



(5.1) 



where, Xc = ((^^(t), ^^(0))) = -^e(^)([nd(^), nd(0)]_) 
denotes the charge susceptibility. Figure [5] shows the 
charge susceptibility Xc as a function of the e-ph coupling 
X/uJo- We see that the charge susceptibility presents a 
wide peak for A = (almost flat charge susceptibility 
spectrum). Then, the broad peak decreases and moves 
towards higher frequencies whereas a low energy peak 
emerges in the spectrum as long as X/ujq grows. The low 
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FIG. 3. (Color online) Electronic spectral density Ad{ijj) for 
different values of X/ojo and % — 0. Parameters are Ed ~ 
-0.05, U = 0.1, 2r = 0.016, = 0, and kgT = 0. 
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FIG. 4. (Color online) Electronic spectral density Ad(Lu) for 
X/iJo = 0.4 and different values of %. This corresponds to 
the weak e-ph coupling regime. Parameters are Ed ~ —0.05, 
U = 0.1, 2r = 0.016, (p = 0, and fcsT = 0. 



energy peak corresponds to charge Kondo correlations 
that have a different energy scale {Tk) than the spin 
Kondo correlations (corresponding to the broad peak). 
This result demonstrates that the e-ph coupling induces 
Kondo correlations in the charge sector reflected in the 
appearance of a low energy peak in Xc- 

In the following we analyze the local density of states 
of the molecule for = both in the weak and strong 
e-ph couphng. We start by showing in Fig. [3] the DOS 
in the absence of the direct channel, i.e., W = for the 
symmetric case £d = —U/2. Here, the main effect of the 
e-ph interaction is to broad the Kondo resonance when 
X/ujQ < 1 and to shrink it when X/loq > 1 (see for example 
the curve corresponding to X/ujq = 1.2 which corresponds 
to the strong e-ph coupling regime). In this case, Tk be- 
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FIG. 5. (Color online) Electronic spectral density Ad{iij) for 
A/cjo ~ 1.2 and different values of %■ This corresponds to 
the strong e-ph coupling regime. Parameters are Sd = —0.05, 
U = 0.1, 2T = 0.016, = 0, and ksT = 0. 



FIG. 6. (Color online) Linear conductance Q as a. function of a 
gate voltage Ed for different values of background transmission 
Tb at fcflT = 0. Parameters are U = 0.1, 2r = 0.016, A = 0, 
and (f — 0. 



comes very small as expected from Eq. (|2.32l) . When 
the direct channel is connected, we see in Fig. U that the 
Kondo resonance is almost unaffected in the weak e-ph 
coupling regime {X/loq < 1). However, there is a remark- 
able feature, the DOS becomes asymmetric. The asym- 
metry in the DOS indicates clearly the electron-hole (e-h) 
symmetry breaking. This lack of e-h symmetry has dra- 
matic consequences in the strong e-ph coupling regime 
where the Kondo effect is due to charge fluctuations. 
Fig. [S] shows the behavior of the DOS when 7& 7^ 0. For 
Tb = we see that the charge Kondo resonance is pinned 
at Ep due to the particle-hole symmetry (remember that 
Sd = —U /2). As long as Tb is turn on, the e-h symmetry is 
broken and the charge Kondo effect is rapidly destroyed. 
All these features are displayed in the linear conductance. 
Here, we have analyzed the linear conductance for differ- 
ent e-ph couplings, different background transmission Tb 
and Aharonov-Bohm phases. Fig. [B] displays the zero 
temperature linear conductance C/ as a function of the 
level position Sd in the absence of the e-ph coupling, i.e., 
A = 0. Different curves correspond to different values of 
background transmission %■ For % — the conductance 
reaches the unitary limit, i.e., Q = 2e'^/h in the Kondo 
regime {—U + T < Ed < — F) and then it drops rapidly 
when charge fluctuations are turned on, i.e., whenever 
-T <ed<T and -CZ-F < Ed < -U + T. Finally, in the 
empty orbital regime the molecule charge rid ~ 0,2 lead- 
ing to a vanishing conductance. This scenario is totally 
modified when the direct path is connected. In such a 
case, interference effects between the resonant and non- 
resonant paths give rise to asymmetric line-shape for the 
conductance and eventually a zero conductance region 
for 7& = 1 in the Kondo regime, due to the destructive 
interference. There is however a tiny feature in the lin- 
ear conductance not noticed so far due to the lack of e-h 
symmetry. The Kondo plateau (for Tb = 1) is not totally 



symmetric with respect to Sd/U = —0.5 as it should be 
in the case of e-h symmetry. This is because the direct 
path breaks this symmetry. However, this has a negligi- 
ble impact when the Kondo correlations are generated by 
spin-fluctuations as we have shown. However, the lack of 
e-h symmetry has stronger consequences on the conduc- 
tance in the strong e-ph coupling corresponding to the 
charge Kondo effect. 

Now, let us discuss the conductance behavior with Tb 
in terms of resonant and non-resonant phase shifts. In 
the Kondo regime (A = 0) we can safely assume that the 
resonant phase shift 5 res = 7r/2. Taking this fact into 
account Eq. (|4.12p is greatly simplified 

T{EF)[5res = ^/2] =a[l- rfcCos2((^)] , (5.2) 

Therefore the linear conductance decreases as long as we 
open the direct path by increasing Tb. It is important to 
realize that Kondo correlations are not destroyed in spite 
of the fact that linear conductance is suppressed. The re- 
duction of the conductance is rather due to interference 
between the lower and the upper paths. Notice, that here 
the conductance reduction is maintained for any value of 
Ed in the Kondo regime {—U+T < Ed < — F), whereas for 
a noninteracting quantum dot this is only true for the res- 
onant condition when the resonant level matches with the 
Fermi energy in which the phase is Sres = n/2. Finally, 
away from the Kondo region Sres 0, tt, and then T{Ep) 
is equivalent to the non-resonant transmission probabil- 
ity Tb- In general for arbitrary 5res and (/3 = 0, we obtain 

T{Ef) = {Tb + aTZb) sir,^{5res - ^„o„) ■ (5.3) 

This expression leads to a conductance peak whenever 
5res — <5non±7''/2 and a vanishing conductance for 5res = 
5non- For Tfc = 0, in the Kondo region, the conductance 
reaches its unitary value {5res = 7''/2 and (5„o„ = 0) and 
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FIG. 7. (Color online) Linear conductance 5 as a function of a 
gate voltage Ed for different values of background transmission 
% at ksT = 0. Parameters are U = 0.1, 2r = 0.016, A/uo = 
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FIG. 8. (Color online) Linear conductance Q versus gate volt- 
age Ed for different values of % at fcsT = in the strong 
e-ph coupling regime, X/too = 1.2. Parameters are U = 0.1, 
2r = 0.016, ujo = 0.05, and ip ^ 0. 



vanishes outside {Sres = Snon = 0). Now we turn on 
Tb and the linear conductance is gradually suppressed in 
the Kendo regime when Tb is augmented [Sres — n/2 and 
cot(5„on) = ^/oRb/Tb, i.e., 6non -> for Tb ^ and 
5non 7r/2 for 7b — > 1]. Outside the Kondo regime, Q 
increases {Sres = and (5„on 7^ 0) . EventuaUy, for Tb = 1 , 
Q vanishes in the Kondo regime {5non,res = ti'/S). Away 
from this regime, the conductance presents a peak since 

Snon = Vr/2 and 6res = 0. 

This physical scenario is modified when A 7^ as 
shown in Fig. [T] The main difference with the A = 
case (cf. Fig. [B]) is that the conductance width decreases 
due to the presence of phonons. This conductance width 
is roughly given by the e-e repulsion, since U is effec- 



tively renormalized by the electron- phonon interaction 
(t^cff — U — 2X^ /uJo) the Kondo plateau region is reduced 
as long as X/ujq increases. A more striking feature ob- 
served in the linear conductance is the evident lack of 
e-h symmetry when 7b 7^ 0. In Fig. |B] when 7b = 1 
the conductance curve is not symmetric with respect to 
£d/U = —0.5 reflecting the lack of e-h symmetry. The ef- 
fect of the e-h symmetry breaking is more pronounced as 
long as A/wo grows and the system enters in the charge 
Kondo effect originated by the bipolaronic attraction. In 
this regime, the lack of e-h symmetry for Sd/U — —0.5 
results in a complete suppression of the charge Kondo 
effect (A/wq > 1) with strong consequences in the hnear 
conductance. The linear conductance for the bipolaron 
regime (strong e-ph coupling limit) is plotted in Fig.[S] In 
the absence of the direct path the conductance consists 
of a very narrow peak of unitary conductance around 
£d/U — —0.5. The width of such a peak reflects the fact 
that the charge Kondo effect occurs only for A w Tk 
(with Tk being the charge Kondo temperature). Any 
deviation from the e-h symmetry point makes the con- 
ductance drops very rapidly to zero. To turn on 7b has 
two different consequences. First, the direct path breaks 
the e-h symmetry and therefore the charge Kondo effect. 
However, by tuning properly Sd/U the Kondo resonance 
can be restored. Second, once the Kondo resonance is re- 
stored at certain value of £d/U , the interference between 
the Kondo resonance and the direct path can occur. 
These two effects are nicely illustrated in Figs. [8] and [H 
Let us consider the linear conductance for 7b = 0.33 (dot- 
dashed curve in Fig. [5]). Here, Q shows a Fano-like shape 
where the maximal and minimal conductance take place 
for a value of Sd/U slightly above the symmetric case 
(edjU = —0.5). Around this value the DOS displays a 
charge Kondo resonance (not shown here). Even more 
evident is the effect of the direct path when 7b = 1- As 
shown in Fig. [8] the complete suppression of the linear 
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FIG. 12. (Color online) Electronic spectral density Ad{ij) 
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FIG. 11. (Color online) AB oscillations of the linear conduc- 
tance C/atfc_BT = 0asa function of ip for % = 1. Parameters 
are Sd = -0.05, U = 0.1, 2r = 0.016, and ujq = 0.05. 



conductance occurs for Sd/U = —0.43 since the Kondo 
restoration takes place at this gate value. The restora- 
tion of the Kondo resonance is showed in Fig. [Q] Here, 
the DOS displays a Kondo resonance at Ep. Neverthe- 
less, the restoration of the Kondo effect does not imply 
the restoration of the e-h symmetry since the DOS for 
£d/U = —0.43 is asymmetric (see in Fig. [9]). 

Finally we comment the results for the linear conduc- 
tance as a function of the AB phase for the two regimes. 
In Fig. [TU] and Fig. [TTl we display the linear conductance 
Q as a function of the AB phase (p. We have considered 
only the symmetric case, ej, = —U/2. From Eq. (j4.12p we 
note that Gi'p) = G{—'p)- Thus, we plot only Lp e [0,7r]. 
As we mentioned the direct path breaks the e-h symme- 
try. This is true for any value of p except for the special 
case of ip = 7r/2. For ip — tt/2 the charge Kondo state is 



not destroyed by the direct channel and the linear con- 
ductance is given by means of Eq. (|5.2p . This is clearer if 
we use the parity basis where the odd and even channels 
are equally coupled to the molecule and the e-h sym- 
metry remains unaltered independently of Tb. To illus- 
trate this we have plotted in Fig. [T^] the spectral density 
at Lp = 7r/2 for various values of the background trans- 
mission Tfe. Having this in mind we can easily explain 
the results plotted in Fig. [101 In the absence of e-ph 
coupling, the linear conductance for the Kondo regime 
follows Eq. (I5.2p . For the weak e-ph coupling regime 
(A/wo = 0.4, and X/ujq = 0.8) Eq. is still vahd to 

describe Q{p>) since for = —0.5 the system is found 
in the purely spin Kondo regime, only U — > f/e// and 
£(i ^d,eff- However, for X/ujq — 1.2 the system is in 
the strong e-ph coupling regime and the Kondo effect is 
partially suppressed for (p ^ n/^ due to the e-h sym- 
metry breaking. For ip k it/2 the conductance shows a 
peak and only for = 7r/2 the conductance is unitary 
due to the presence of the charge Kondo resonance and 
its constructive interference with the non-resonant path. 
For Tb = 0.33 and (/? « 0, tt the Kondo resonance is not 
totally quenched (see Fig. [5] for 1^9 = 0) and therefore 
the conductance is not only reduced to the background 
transmission but it has some contributions coming from 
interference between the Kondo resonant path and the 
direct path. 

Fig. [TT] shows Q versus Lp when the direct path is com- 
pletely opened, i.e., % = I- For A/wq = 0,0.4,0.8 all 
these curves follows Eq. (|5.2p that results from the in- 
terference between the spin Kondo resonance and the di- 
rect channel. However for X/luq = 1.2 the conductance 
shows only contribution of the background transmission 
probability and a tiny contribution attributed to the in- 
terference between a strong suppressed Kondo effect and 
the direct path that makes not unitary the conductance 
for « 7r,0. At (p = Tr/2 the conductance reaches the 
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unitary value since for this special phase value the charge 
Kondo effect is not destroyed and we have constructive 
interference. 



VI. CONCLUSIONS 

We have analyzed the linear transport properties of an 
AB interferometer with a molecule inserted in one of its 
arms. The vibronic internal degrees of freedom of the 
molecule are described by means of phonons modes cou- 
pled to the molecular occupation. Depending on the e-ph 
coupling we distinguish two different scenarios, the weak 
e-ph coupling limit where the electronic transport occurs 
as single electron tunneling events and the strong e-ph in- 
teraction which leads to bipolaronic attraction between 
electrons. Here the tunneling through the molecule takes 
place in pair of electrons. We investigate the low tem- 
perature limit where Kondo correlations emerge when 
the molecule is strongly coupled to the leads. Depend- 
ing on the e-ph coupling the Kondo correlations have a 
different origin. For the weak e-ph coupling the Kondo 
effect corresponds to the spin Kondo effect, while for a 
strong e-ph interaction the bipolaron formation leads to 
the charge Kondo effect with fluctuations between zero 
and double charge states (bipolaron state). The effect of 
the direct path breaks the e-h symmetry with dramatic 
consequences in the transport properties. We show that 
the direct path destroys completely the charge Kondo 
effect and suppresses the interference between the res- 
onant and non-resonant paths. By a proper tuning of 
the molecular level position, it is possible to restore the 
Kondo resonance but not the e-h symmetry. In this case, 
when the charge Kondo state is revived, interference ef- 
fects come up. The usual spin Kondo effect is less sen- 
sitive to the e-h symmetry and is therefore less affected 
by the direct path. In this case, the transport proper- 
ties in the weak e-ph coupling regime are dominated by 
the interferences between the spin Kondo effect and the 
direct path. In order to obtain these results, we have 
chosen values for the parameters which are in the range 
of typical experiments in molecular transport and we are 
therefore confident that our predictions can be tested in 
the near future. 



ACKNOWLEDGMENTS 

We acknowledge D. Sanchez and R. Zitko for their 
fruitful discussions. PS would like to acknowledge O. 
Entin-Wohlman and A. Aharony for an earlier collabo- 
ration on related issues. R.L. and J.S.L acknowledge fi- 
nancial support from grant FIS2008-00781 and the Con- 
selleria d'Innovacio , Interior i Justicia (before Consel- 
leria d'Economia, Hisenda i Innovacio), Govern de les 
Illes Baleares. G.P acknowledges finnantial support from 
grant: MEG: MAT2008-02626. 



Appendix A: Mapping to a Kondo Hamiltonian 

In this appendix we give the details for the derivation 
of effective Hamiltonian using the projection method. 
This consists of a second order Rayleigh-Schrodinger 
perturbative calculation on the tunneling Hamiltonian, 
We analyze the two cases, namely, the weak e-ph 
coupling and the strong e-ph coupling regime. 

Let us take ipn as the component in which the occupa- 
tion number of the molecule is rid, then we can write the 
following matrix derived from the Schrodinger equation 



(Al) 



PnHPn' ■ Here, P„ are the local projectors 



E — Hqo 


—Hoi 


— ^02 




ipo 






E-Hii 


— ^12 










—'H21 


E-n22 




1p2 





where Hnn' 
defined as 



where 



Po = |0)(0| =(l-nt)(l-^i): 
^t = lt)(tl=«t(l-"i): 

n = u)ai = "i(i-"t): 
P2 = \n){n\ ="T"i- 



Po+Pt + n+^2 - 1. 



(A2) 



(A3) 



In this fashion the effective Hamiltonian will be given by 
(we only write the off-diagonal contributions) 



n 



eff 



1 



E 



1,0 



E, 



-Hn 



0,m 



El n — E- 



2,m 



-H21, (A4) 



where the energies -Ei.o, £'0,m, £'2,m are the eigenener- 
gies of the isolated molecule coupled to the local phonon 
mode. 



1. Weak electron-phonon coupling limit 

The eigenenergies of the isolated molecule and the lo- 
cal projectors to perform the second-order perturbative 
calculation are given respectively by 

Ei,o = Sd , 

Eo,m = 1- mujQ , 

Wo 
X2 

E2.n1 = V2ed + U + mujQ , 

Wo 

Hoi = \/2 Vtimp-^cl^J^il-Uda), 

i,k,cr.7n 

H2i = V2 J2 V;{m\U+\0)sgn{a)dindaCika , 
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with a = —<j. The tunnehng amplitudes for the even 
and odd channels are respectively, K = cos((/?/2)F and 
Vo = sin((p/2)V. Using the above expressions the effec- 
tive Hamiltonian is 



p w ^ ^ X! X! X! 

^1,0 — -C'O.m 



y„*yff|(m|^Y-|0)p 

\ 2 



(1 - naa)dlcapaclq„,d„' (1 - "-dff' ) , 

' a,p,fT p^q,a' m ujq " ^ 

sgn(cr)sgn(cr')cJjp^?i<j<TdCTrfI'"d<T'C;3q<T' , (A6) 

Summing over the spin indices all these contributions 
yields 

a.,/3,p.q 

- 51 J2apndidfd\ndiclp^ci3g^ , a =t, cr' =t, 

- X! '^Oa/3(l - ndl)d^-fdl{l - ndf)Capfcjjg^^ 

- ^ J2a0ndidfdlnd^clp^ci3qi , a =t, o"' =1, 

- ^ J2a0nd^did\ndiclp^ci3g^ , a =1, a' =t, 

a,/3,p,q 
a,/3,p,g 

where we have defined the following quantities as the 
exchanged couplings 



Using these results Eq. (jA4p reads 

"^eff — ^ JQal3 ?^dt ( 1 ~ '^-d^) (1 ~ <5a,^(5p,q)c|3^^Cap^ 

+ d\dlCp^^Capf + d'^^dfc'p^^Capl 
+ «rf4,(l - «rft)(l ~ Sa,pSp^q)c^p^^Capi 



0.^/3 ^p^q 



2V*Vp\{m\U-\0)\' 
- — Sd — rncuQ 



(A7) 



(A9) 



Since the possibility of p = g is extremely small, we can 
safely drop the Kronecker delta terms. By employing the 
following identities 



5"+ = d\di , S- = dldf , 
4^. jt 



(AlOa) 



(AlOb) 

1 = "-dt + ridi = rid-fil ~ Udi) + Udiil - Ud-f) . 

(AlOc) 



Equation (|A9[) can be then written as 



+ ■ clp;C/3,^ + S ■ clp^C/Sqi 



Cpqa 



(All) 



where J"'' = ( Joa/3 + J2a/3) and /ST"'' = ( Joa/3 - J2ap) 



2. strong electron-phonon coupling: the 
bipolaronic scenario 



2V^V^\{m\U+\0)\^ 



J, 



2a/3 



(A8) 



Under the particle-hole transformation, Eq. I\2.7\i be- 
comes 



£=l/2,fe,(T 



(A12a) 

+ E [sgn{a)We''^cl^,^cika + h.c] , (A12b) 
■Hji/ = E £dad\d„ + Uud^ridi , 



{^'^kaclk^d^ + h.C. 

£=l/2,fe,(T 

Ue-ph = A(a + a^) {nd\ - Udi) , 



(Al2c) 

(A12d) 

(Al2e) 
(Al2f) 
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where Sda = U/2 + sgn(cr)A/2, U = -U and A = 
2ed + U. Assumming the symmetry condition for the 
conduction band Eka- = ^^ka^ tunnehng ampUtudes 
become 



(A13) 
(A14) 



Now we perform the Lang-Firsov transformationj^^iS 
The molecular part in the presence of phonons {T-Lm + 
He-ph + 'Hph) in Eq. (|A12I) reads 



(A15) 



+ f [/ + — — j ridfUdi + ujQtt^a , 
and the tunneling part "Ht is 

i=l/2,k,a- 

Next we perform a canonical transformation of the total 
Hamihonian {H O'^HO) using the following unitary 
operator 



O = e~' 



(A17) 



This transformation yields the following eigenstates and 
eigenvalues for the isolated molecule (in the equivalent 
model): 

A2 



|0,m) = |0)|m), = \- muo 



(A18a) 

LUQ 

5,m) =Z^^g"(^)|5)|m), Es,^ = Ida + mcjo, (A18b) 

|2,m) = |2)|m), E^j.m = — + Edo- + ed<f + + mwo, 

(Al8c) 



with 



WsnM = exp 



Wo 



sgn(cr) — (a'l' - a) 



(A19) 



This define the projectors for the strong e-ph coupling 
regime as 

Hoi = V2 ^ i^tt.(m|Z^-^s"(^)|0)4^d„(l-nd<,), 



(A20a) 

H2i = V2 V,laH^^'^^^"^\0)sgn{a)dind^c,ka, 

(A20b) 



i,k,{T,7n 



with cr = -cr, and V^fccr = cos((p/2)Vfccr and Voka- = 
sin((y9/2)VA;cr- With these ingredients the effective Hamil- 



tonian is given by the sum of the following terms 
1 



EE 

a,p,(7 I3,q,a' in 



■n 



01 



A2 



(1 - nds)dlcapaclqcr''^a' ~ '^da') (A21a) 



^1 



-n2 



2.m 



= 2E EE- 



sgn{a)sgn{a')clp„ndadcr3l,nda'Cpqcr' ■ (A21b) 

The energies Ea-^ of the initial and final states can be 
different. We eliminate this ambiguity by choosing the 
average between both energy states. Finally, a similar 
calculation than for the weak regime, produces the fol- 
lowing exchange couplings 



,a/3 _ _8 ^ VaV0\{m\U-\O}\^ 



u 



„/j _ ^ V- VaV0{O\U-\m){m\U-\O) 



_2Ai 



2mLJo 
U 



(A22a) 
(A22b) 



U ^ 



2A3 



1 



Now we can perform an expansion in the large A/wq limit 

2 



J,? 



V 



= -(A/i^o)' 



_2Ai 



1 



2A^ -I I 2aJo 
wo(7 V 



2A2'' 
2A2 ' 



jal3 



u 



1 / A 



2! V wo 



-(A/wo) 



2A2 



1 



2A^ -I I 2uio 



2A^' 

woC/ 



-(AM)^ 



-2(A/t^o)' 



i_ I AV + 1 f A 

Loo) 2! Icjo 



and finally we note that is exponentially suppressed 
in comparison with Jy giving a very small Kondo temper- 
ature in the strong e-ph coupling regime [see Eq. (I2.30p ] . 



Appendix B: Electronic and phonon Green functions 



In this appendix we calculate the electronic and 
phonon Green functions in order to calculate the den- 
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sity of states of electrons and the linear conductance. 
To this end we employ the equation of motion (EOM) 
method.— In this manner the equation of motion for the 
correlator < A, B > in the frequency domain (z) is given 

by 

z{{A,B)), + {{[H,A]^,B)) = {[A,B]i), (Bl) 

with = + if both A and B are ferniionic operators, 
while ( — — otherwise. Notice that we have defined as 

/oo 
dte'^'{{A,B))t, (B2a) 
-OO 

{{A,B))t = -te{t){[A{t),B]^). (B2b) 

Here, 8 denotes the step function, (•) is the thermo- 
dynamic average. Next, we obtain the electron Green 
function in the presence of the e-ph interaction. We 
derived the molecule Green function and also the lead- 
molecule Green function. For A = and A — Caka 
with B = d\, the equation of motion for the molecule 
and lead-molecule Green function are respectively 

[z ~ ed)gda,d'y{z) - 1 + U{{d,dUs,4)), (B3) 
+X{{d,{a + at), 4)), + ^, F„*((c„,^, 4)), , 

and 

+ p^e-gn(a)^^^^^^^^^^t))^^ (B4) 
9 

where sgn(a) = — (+) for a = L{R) and a — R{L). Here, 
Eq. (|B4|) can be rewritten as 



this we write them using the spectral representation 

{{dj^d-,,dl)), 
= -:^ ^{n\dl\m){m\d^dlds\n) {e' 

1 



V 



Z - {En + Em) 



- in5{z - En + E,n) 



, (B8) 



and 



((d,(a + at),4))^ 

= ^^(r^|4l'7^)(TO|(i^(a^-at)|n) 



V 



Z - {En + Em) 



iTTS{z — En + E„ 



, (B9) 



where Z is the partition function and V denotes the prin- 
cipal value in the Cauchy sense. 

The phonon Green function is calculated in the same 
manner as the electronic Green function by using the 
EOM method. Now A = a and £? = and then the 
phonon Green's function reads 



(z - wo)((a, a^))z = 1 + A((nrf, a^)), 



(BIO) 



By dividing both sides of Eq. (|B10p by {{a,a^))z we get 
a Dyson type equation for the phonon propagator 



V-\z)=V,\z)-\ 



(Bll) 



p 

(B5a) 

q 

(B5b) 



where I'(z) = {{a,a^))z, Vq^{z) = z — ujo and Ud — 
Y.a'^da - 1- The propagators {{nd,a^))z and {{a,a^))z 
are obtained from the NRG spectrum. Using again the 
spectral representation we have: 



((«'«^)>^--4Eih«^ 



-I3E„ 



where ((c„i„4)), = Ql^i.^dai^)- Taking Eqs. dUil, 
(|B5bp . Eq. (|B3l) . and Eq. ((Bi)) the molecule Green func- 
tion reads 



Z — {En — Em) 



- inS {z - {En - Em)) 



) 



(B12) 



and 



- X 



{{d,dUa,di))z 
Qda,dc7{z) 

((d,(a + at),4)). 



(B6) 



{{nd,a^))z = -\ E("l«V)("^l^<i|«) (e-^""" - e"^^'") 



Qda,da{z) 



Z — {En — Em) 



- inS {z - {En - Em)) 



(B13) 



where 



^dJ,d.i^) = z^ed-^[-2ir-27T'pl\V\'Wcos{ip)] . 

(B7) 

Finally, {{dcd%ds,dl)) z and {{da{a + a^),dl)) z are de- 
termined from the NRG energy spectrum. In order to do 



This method of calculating the phonon propagator 
through the EOM technique [see Eq. (jBllI) ] is much more 
accurate than the direct calculation of this from the NRG 
spectrum [see Eq. (|B12[) ] as it has been shown by Jeon 
et at, in Ref. M- 
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